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Abstract 

We study the non-equiliribium dynamics of atom-molecule Bose gases in a double-well potential. 
In this system, the internal atom-molecule tunneling has significant influence on the dynamics. 
We investigate the periodicity of dynamics by studying the level statistics of the quantum system. 
We find that chaotic energy eigenstates arise from the competition between the interwell and 
the atom-molecule internal tunnelings. Furthermore, we show that the physical quantities relax 
to the micro canonical averages in the full-quantum dynamics when the system is chaotic. This 
thermalization is caused by the verification of eigenstate thermalization hypothesis (ETH). We 
show numerically that the onset of ETH occurs simultaneously with that of chaos. In addition, we 
show that the energy eigenstates become to be exponentially localized states simultaneously with 
the onset of chaos. 
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I. INTRODUCTION 



Recently, the experimental studies of quantum signature of chaos have been achieved 



le dynamical localization, which is the 



6|, has been observed. The quantum 



in ultracold atomic systems [miM- For instance, t 
quantum counterpart of chaos in classical systems 
signature of chaos has become one of the important subjects in cold atom physics. 

In this paper, we study the non-equiliribium quantum dynamics of atom-molecule Bose 
gases in a double-well potential. Recently, the tunneling dynamics between left and right 
wells is observed in Bose- Einstein condensates trapped in a double- well potential [3]. This 
experimental achievement has triggered much interesting research Ulioi. Though the spa- 
tial coherence of Bose gases is focused in these studies, tunneling effects occur not only 
between spatially separated states but also between internal degrees of freedom. In particu- 
lar, the effects of internal tunneling between atomic and molecular states have been discussed 



theoretically 



llNl3|. Furthermore, in an atom-molecule mixture in an optical lattice, the 
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Iq . Based on these stud- 



atom-molecule internal tunneling dynamics has been observed 
ies, we study atom-molecule Bose gases in a double-well potential with a particular focus on 
the effect of atom-molecule internal tunneling. As a result, we show that the atom- molecule 
internal tunneling induces chaos and thermalization dynamics. 

In general full-quantum systems, the chaotic dynamics can not occur because the 
Schrodinger equation is the linear equation (see Ref. and references there in). The 



quantum signature of the semiclassical non-periodicity is an important problem 18 



. In the 



preceding studies [l9|, |20|, it is conjectured that the quantum signature of chaos appears 



in the level spacing distribution. Based on these considerations, we investigate the statis- 
tical property of energy spectra. As a result, we conclude that the atom- molecule internal 
tunneling induces the chaotic energy eigenstates. 

Furthermore, we discuss the relation between chaotic energy eigenstates and thermaliza- 
tion. Recently, the experimental challenge for the foundation of statistical mechanics has 
become possible. By using quasi-one-dimensional Bose gases, the non-equilibrium dynamics 



of isolated quantum many-body systems is investigated experimentally 
to 250 ^''Rb atoms do not thermalize 



21 



22j. From 40 



2l|. It is showed that the multiplicity of degrees of 



freedom induces thermalization generally 



23j. In the experiment 2l|, it seems that the 



number of degrees of freedom is not large enough for thermalization. 
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However some numerical studies exhibit that the non-integrabihty of systems induced 



thermahzation even in a few body systems 



24H30| ■ In these studies, the integrabihty of 



the Bose-Hubbard model is estimated by investigating the level spacing statistics. When 
the system is non-integrable, the semi-classical dynamics is chaotic, otherwise the system 
exhibits regular dynamics. The preceding studies investigated the quantum dynamics in the 
non-integrable region. The numerical results indicate that when the semi-classical dynam- 
ics is chaotic, thermahzation occurs in the corresponding quantum dynamics. Otherwise 
thermahzation does not occur. The integrabihty of systems plays an important role in the 
thermahzation mechanism of isolated-quantum systems. 

In addition, the system size is also important for thermahzation. The preceding studies 
focused on single-component Bose gases in a optical lattice. These studies showed the 
appearance of non-integrability by investigating the Bose gasses confined in about 20 sites 
and 8~12 sites |30(] optical lattices. In particular, KoUath et al. 30| showed that the 



non-integrability is increased by increasing the system size. From this result, the system 
size should be large to some extent in order for thermahzation to occur. Furthermore, in the 
preceding studies, in the region where the interwell tunneling and interparticle interaction 
compete, the system exhibits thermahzation. 

In contrast to these studies, we show that thermahzation can occur in the system with 
only two sites. Our subject is atom-molecule Bose gases in a double-well potential. In 
this system, the internal atom-molecule tunneling induces chaos. In addition, we confirm 
that thermahzation occurs in the corresponding quantum system. The large system size is 
not always needed for thermahzation. Furthermore, this thermahzation process dose not 
need interparticle interactions. The competition between the interwell tunneling and the 
atom-molecule internal tunneling induces the thermahzation dynamics. 

We also find that the energy eigenstates become to be exponentially localized. This is 
quite similar to Anderson localization. In the previous studies, it is pointed out the relation 
between quantum chaos and the Anderson localization in Kicked rotor 5[. Furthermore, 
this is confirmed clearly in the experiments 2|, |3[. However, this relation is not confirmed in 
general cases. In this paper, we investigate whether the exponential localization occurs in 
our system, and find it does indeed. 
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II. MODEL 



A. Hamiltonian 

The second-quantized Hamiltonian for Bose atoms and molecules can be written as 



where \E'a and represent field operators for Bose atoms and molecules respectively, A is the 
internal tunneling strength between atomic and molecular states, S is the energy difference 
between atoms and molecules, and Kxtlr) is a double- well potential, rria is the atomic 
mass, rrih = 2ma is the molecular mass. The interatomic interaction can be approximated 
in terms of the s-wave scattering lengths as Qa = Airh^as/nia. The intermolecule, and the 
atom-molecule interactions are neglected in our treatment. Furthermore, we introduce the 
four-mode approximation. In this approximation, we concentrate on only lowest-energy 
atomic and molecular modes in each well, and ignore the effect of the particles occupying 
other modes. From this point of view, field operators can be approximated as — ^abO-L + 
^aRCiR, ^6 ~ ^hhbL + ^bRhR, where ^aR i^bL, ^bR) are the wavefunctions of the atomic 
(molecular) lowest-energy modes in the left well and the right well respectively. ai^aRipL, &r) 
are annihilation operators for the atomic (molecular) lowest-energy modes in the left well 
and the right well respectively. Applying these approximations to Eq. ([1]), we obtain the 
quantum four- mode Hamiltonian (four- mode model). 



where the parameters are defined in Appendix. In order to study the non-equilibrium 
dynamics, we obtain all the eigenenergies and eigenvectors by performing the full exact 
diagonalization of the Hamiltonian ([2]). 



^ = E hA ^^^^^ ■ + Kxt(r)^l^. + f / ^r^i^l^-^ 




a 



(1) 



H = -Ja{a{aR + aiai) - Jb{bibR + + A{b[bL + 6^6^) 

I ^» (' t t I t t N 

+-^{(^L(^L(^LaL + aRdRaRaR) 




(2) 
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In what follows, we write the ith energy eigenstate as 

= Yl (^"i' ^-^^ ^bL, NbR) \NaL, NaR, NbL, iVfcij), (3) 

where NaL{bL) and NbLi^hR) are the atomic (molecular) particle numbers in the left and right 
wells. \NaL, NaR, NbL, Nijr) is the eigenstate of the operators o'^cil, o^j^aR, h^hi and 
The total particle number = NaL + NaR, + 2 (A^'^^ + A'';,^) is always conserved. In this 
paper, we calculate the summation YliNaL NaR n^l n^r wilder the condition of the conservation 
of the total particle number A^. We write the expectation value of the ith energy eigenstate 
as {A)i = A is any observable. 

We set the parameters in the Hamiltonian. We set the molecular-tunneling strength as 
Jb/Ja = |- Next, we consider the atom-molecule energy difference A. This is controllable 
parameter in the experiment by means of Feshbach resonances. In this study, we set this 
value as A/J^ = —1. We set Ua = except Sec. |VT1 In only Sec. IVTI we set Ua and 
investigate the influence of the interparticle interaction on our results. 



B. Time-evolution of particle numbers 

The initial state is expanded as (t = 0)) = J2a Cal'^a), where |\E'a) represents the ath 
energy eigenstate. The many-body wavefunction evolves as (t)) = e~^^*^'^\'^ (t = 0)) = 
Cae~*^"*|\E'Q,), where Ea is the a's eigenenergies. The quantum-mechanial expectation 
of any observable A evolves as 

(i) = (vl> (t) |i|vl> (t)) = ^C:C;3e-<^"-^^)*A«^, (4) 
where Aai3 = {'^ alA]"^ l^) ■ The long-time average of (A) is 

'{A) = J2\Ca\'Aaa. (5) 
a 

If the system relaxes to equilibration, this is the relaxed value of the physical quantities. 
In this paper, we investigate whether the physical quantities relax to the microcanonical 
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ensemble average. The microcanonical ensemble average is given as 

i'=i-AN,i+AN 

where the ith energy eigenstate has the closest energy eigenvalue to the mean energy of the 
initial state. We take the 2 x AA^ energy eigenstates around the i'th energy eigenstate, and 
calculate the microcanonical average using the principle of equal weight. Unless AA^ is not 
too large, any other values of AA^ does not change our results as explained in Sec. IIVBI 

III. LEVEL STATISTICS 

In this section, we show that the atom-molecule internal tunneling induces the chaotic 
energy eigenstates. First, we introduce the indicator for the quantum signature reflecting 
chaos in corresponding semi-classical systems. The quantum signature of chaos appears in 
level spacing distributions. 

Here we explain how the author calculates level spacing distribution. First we perform 
the exact diagonalization of the Hamiltonian matrix (|2]) and obtain the eigenenergies Ei 
{i = 1,...,Nt)- Nt is the number of the eigenstates. Ei are arranged in ascending order. 
Second we calculate the adjacent level spacings as Si = Ei^i — Ei {i = 1, ...,Nt — 1). We 
define the average of the adjacent level spacings as S* = J2f=]~^ ^-il i^T — !)• Next we divide 
the {Enj, —El) by 0.16 x S. Then we count the number of eigenstates included in the intervals 
A^i = [(z - 1) X 0.165, z X 0.16S']. Finally we write the histo gram. The horizontal axis is 
normalized by S", and the vertical axis is normalized so as to / dSP {S) = 1. We note that 
we can change 0.16. Our results are not so sensitive to this value. 

In a quantum mechanical system, the dynamics is always periodic and chaos does not 
occur, because the time-evoluiton equation is linear. However, the breakdown of periodicity 
in the classical counterpart has influence on the level spacing distribution of the quantum 
system, because the level spacing is related with the dynamical nature of physical quanti- 
ties as {A) = 'Yl,a i3^oi^P^~^^^°'~^'^'^^ ^oip- Intuitively, level repulsions correspond to chaos. 
When chaos occurs, many periods are included in the dynamics in the semiclassical case. 
Therefore, the transition between energy eigenstates occur in the quantum counterpart, and 
level repulsions are induced. In particular it has been conjectured that, when chaos occurs 
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FIG. 1: (color online) Level spacing distri- 
butions at \fNg/Ja = 0.5. = 40. Blue 
and violet bars represent odd and even parity 
energy spectra. Red and green lines are Pois- 
son and Wigner distributions. The lower and 
higher 20 % levels are not included. 



FIG. 2: (color online) Level spacing distri- 
butions at y/Ng/Ja = 1.5 and = 40. Blue 
and violet bars represent odd and even par- 
ity energy spectra. Red and green lines are 
Poisson and Wigner distributions. Lower 20 
% and higher 20 % levels are omitted. 



in the classical counterpart, the level spacing distribution in the quantum system shows the 
universal features, which are related to the random matrix theory. 19 1 

From the random matrix theory, the level-spacing distribution of a spinless and time- 
revaersal Hamiltonian exhibits the Wigner-Dyson distribution 



^H^(5) = ^exp(-7r5V4), 



(7) 



where S is the adjacent level spacing normalized by the mean-level spacing. Pw [S] takes the 
maximum value at 7^ 0. This indicates that level repulsions occur. When the dynamics 
is chaotic in classical systems, it is conjectured that the level-spacing distribution of the 
quantum counterpart is the Wigner-Dyson distribution 19|. In contrast to chaotic cases, 
when the system is integrable, the level-spacing distribution is the Poissonian distribution 



3l| 



Pp{S) 



(8) 



If the level-spacing distribution follows the Poissonian distribution, the semi-classical dy- 
namics is regular. In contrast to chaotic cases, Pp (S) has the maximum value at S* = 0. 
This indicates that level clustering occurs. 

Before the analysis, we need to classify the energy spectra according to the symmetry 
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of the system. The symmetry has significant influence on level repulsions [32|. Since the 
double-well potential has the left-right symmetry, the energy eigenstates can be classified 
into odd and even parity states. By interchanging the left and right well particle numbers of 
the eigenfunctions as NaL ^ NaR and N^^ o ^feR, even parity states are not changed, while 
the sign of odd parity states are changed as — |^)- We should investigate the odd and 

even parity spectra separately. In Figs. [1] and [21 we show the level spacing distribution of 
the even and odd parity spectra. S is normalized by the mean adjacent level spacing S. We 
do not include the lower 20 % and higher 20 % levels because these eigenstates are always 
regular, not chaotic. 

When the atom-molecule internal tunneling is small {s/l^g/Ja = 0.5), the level spacing 
distributions are close to the Poisson distribution. However, the level spacing distribu- 
tions come close to the Wigner distribution in the region where the atom-molecule internal 
tunneling and the interwell tunneling compete {y/Ng/Ja = 1.5). The dynamical property 
is significantly changed by the atom-molecule internal tunneling. In the next section, we 
investigate the time evolutions of the physical quantities in these regions. 

In order to evaluate the periodicity of semiclassical dynamics quantitatively, we define 
the indicator 



V = 



C dS [P (S) - PwiS)] 



J^''dS[PpiS)-P^iS)] 



(9) 



where Pw (S) and Pp (S) intersect at Sq ~ 0.473. When r] = 1, P (S) = Pp (S) and the 
dynamics is periodic. When rj = 0, P (S) = Pw (S) and the dynamics is chaotic. In Fig. 
[31 we investigate rj by increasing the atom-molecule internal tunneling strength g. In this 
figure, we use the even spectra and do not include the lower and higher 20% energy levels, 
because these energy-eigenstates are always regular. From this figure, we conclude that the 
atom-molecule internal tunneling changes the dynamics drastically. Chaos is induced by the 
atom-molecule internal tunneling. 

Here, we note that our results in Fig. [Sjis independent of the interwell molecular tunneling 
Jf). In the Wigner region (y/Ng/Ja = 4), we investigate rj by changing Ji, from to 1. As 
a result, r] is almost insensitive to Jb. Chaos is induced by the internal atom-molecule 
tunneling independently to the molecular interwell tunneling strength. 

Furthermore, in Fig. [3l we found the scaling law where i] depends on the total particle 
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FIG. 3: (color online) The ^-dependence of rj. In this figure, only even spectra are included. The 
lower and higher 20 % levels are not included. 
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FIG. 4: (color online) Level spacing distribu- 
tions at VNg/Ja = 30 and N = 40. Blue and 
violet bars represent odd and even parity en- 
ergy spectra. Red and green lines are Poisson 
and Wigner distributions. 
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FIG. 5: (color online) The (7-dependence of 
T]. N = AO. The lower and higher 20 % levels 
are not included. 



number with the combination yNg/Ja . This indicates that the quantum fluctuation 
is weak as explained below. In the semi-classical limit, the creation- annihilation operators 
are replaced by c-numbers as ciLiR) ^ A/iVaI(^e*^''^(''«) , b^R) ~ A/iVbI(^e*^''-f'(''«) , and we 
obtain the semi-classical Hamiltonian as 



Hci/N 



-2Ja\/XLXR COS {QaR - ^ ah) " '^Jh^/VLVR COS {9bL " ObR) 

-2\fNg [xl^/vlcos {26aL - Obh) + XRy/y^cos {29aR - 9bR)] , 



(10) 
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where the normahzed particle numbers are defined as XL(ii) = NaL{aR)/N and yL{R) = 
NbL{bR)/N. In the classical Hamiltonian f fTOj) . the atom- molecule internal tunneling strength 
g is scaled by y/N. Therefore, when the quantum fiuctuation is weak, the physical quantities 
depend on g via combination \/Ng. Then we conclude that the system behaves classically 
in a certain sense, although we perform the full-quantum analysis by using the full-exact 
diagonalization method. 

Next, we investigate the large g region. In Fig. HI we show the level spacing distributions 
at y/Ng/Ja = 30. The level spacing distributions are close to the Poisson distribution. In 
Fig. we investigate the (^-dependence of t] in the large g region, t] comes close to unity in 
this region. The high atom-molecule internal tunneling recovers the periodicity of dynamics. 

In the region y/Ng/Ja ^ 1 {y/Ng/Ja ^ 1) regions, the system becomes two independent 
interwell (internal) boson Josephson junctions. In a boson Josephson junction, the dynamical 
variables are the relative particle number and the relative phase. The number of degrees of 
freedom is two. 

In order for chaos to occur, the dimension of the dynamical system must be greater than 
two. Since the trajectories in the phase space cannot intersect with each other, significantly 
complex trajectories cannot exist in 2D systems. However, in higher dimensional systems, 
the trajectories in the phase space travel around much more freely than 2D dynamical 
systems. 

In our system, chaos can occur in the region where the internal and interwell tunneling 
compete. We note that this is the necessary condition but not the sufficient condition for 
chaos. As shown in Fig. [3l chaos does not occur at y/Ng/Ja = 0.4, although the interwell 
and internal tunneling strengths compete. 

IV. THERMALIZATION 
A. relaxation dynamics 

In this section, we investigate the relaxation dynamics in both the Wigner region 
[y/Ng/Ja = 4) and Poisson region {y/Ng/Ja = 0.2). In particular, we focus on the dy- 
namics of the atomic (molecular)-particle numbers in the left and right wells. 

The important nature of thermal equilibration is that the final equilibrium state does not 
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FIG. 6: (color online) The time evolution of atomic particle numbers from the easel initial state. 
We investigate the Wigner region {^fNg/Ja = 4). 
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FIG. 7: (color online) The time evolution of molecular particle numbers from the easel initial 
state. We investigate the Wigner region {^fNg/Ja = 4). 



depend on the initial state. To confirm this, we investigate the relaxation from two different 
initial states, which have the same total energy. We call these different initial states easel 
and case2. In the initial state of easel, total particles localize in the left-well-molecular state 
only as |0, 0, N/2, 0). In the initial state of ease2, a half of particles localizes in the left-well- 
molecular state, and another does in the right-well-molecular state as |0, 0, A^/4). In 
these initial states, the total energy is A^A from ([2]). 

In Figs. El [71 [H] and ini we investigate the time-evolution in the chaotic region {\fNg/ Ja = 
4), and compare the results from the different initial states easel and ease2. As clearly 
shown in these figures, the dynamics is irreversible, i.e. the relaxation dynamics occurs. 
Furthermore, the relaxed values are the same in both the easel and ease2. These relaxed 
values coincide with the microcanonical averages calculated by Eq. ([6|). When we calculate 
the microcanonical averages in this section, we set as AA^ = 20. 
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FIG. 8: (color online) The time evolution of atomic particle numbers from the case2 initial state. 
We investigate the Wigner region {^/Ng/Ja = 4). Red and green lines overlap. 




10 

Wh) t 



20 



FIG. 9: (color online) The time evolution of molecular particle numbers from the case2 initial 
state. We investigate the Wigner region {\fNg/Ja = 4). Red and green lines overlap. 
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FIG. 10: (color online) The distribution 
in the energy-eigenstates in the initial states 
easel ". This data is in the chaotic region 
^/Ng/Ja = 4. Eq is the energy of the ground 
state. 
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FIG. 11: (color online) The distribution 

in the energy-eigenstates in the initial states 
easel " and ^^case2 ". The red represents 
easel ", and the green line is "case2 ". This 

data is in the regular region y/Ng/Ja = 0.2. 

Eq is the energy of the ground state. 
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FIG. 12: (color online) The initial state is easel ". This data is in the regular region yNg/ Ja 
0.2. The microcanonical expectation value is (^|,(^)^L(i?))™"° — 6 in this parameter. 



In contrast to these relaxation dynamics, the dynamics are periodic in the regular region 
\/Ng/ Ja = 0.2. The time evolutions of the molecular particle numbers are shown in Fig. [121 
In this figure, the dynamics is periodic, and the relaxation to the microcanonical average 
does not occur. 

In the last of this section, we note that relaxation is induced by the dephasing process as 
(i) = Y.ClCpe<''-'^^yA^p ^ (11) 

a,/3 a 

In the Wigner region, the distribution in the energy-eigenstates are spread out because of the 
level repulsion. As shown in Fig. [TOl the distribution in the energy eigenstates is close to be 
continuous. In this region, the many different period oscillations are superimposed, and the 
dephasing occurs. In contrast to this chaotic cases, in the Poisson region, the distribution 
in the energy eigenstates is discrete as shown in Fig. [TTl The energy levels are clustered 
such as Pp (S) = . In this region, the relaxation can not occur. 



B. Eigenstate thermalization hypothesis 



In this section, we verify the eigenstate thermalization hypothesis (ETH) [2J]. In Figs. 
El [7l [HI and M The relaxation to the same value occurs, although the distributions in the 
energy eigenstates are different in the initial condition easel and case2, as shown in Figs. 
[131 and [HI This is caused by the significant behavior of the energy eigenfunctions as shown 
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in Fig. [151 In this figure, the quantum-mechanical mean values of physical quantities have 
the same value over many different energy eigenstates. This means that Aa,a is independent 
of a in (ITT]) . This behavior of energy eigenfunctions is called the eigenstate thermalization 
hypothesis. ETH explains the thermalization mechanism in Figs. El El 13 and [HI 

Because of the independence of A^^a on the energy eigenstate labeling a, the constant 
value of Aa^a is almost the same as the microcanonical value defined in Eq. ([6]). We represent 
this value as (^A.)'^^'^^" . Therefore the relaxed value is represented as 

(i) ^ Yl \Ca\'A^,^ ^ (A)™ J] = (A)™. (12) 

a a 

As explained above, when the eigenstate thermalization hypothesis is verified, the physical 
quantities relax to the microcanonical mean values, calculated by the principle of equal 
weight. It is clear from Eq. f|T2|) that the relaxed values are independent of the coefficient 
Ca- Thermal equilibration is independent of the initial states. Therefore, the relaxation 
values are almost the same in the easel and ease2, although the distributions in the energy- 
eigenstates are different as shown in Figs. [T31 and (TH In addition, it is clear from Fig. 
[151 that the microcanonical ensemble average of atomic particle numbers defined as ([6]) 
do not depend on AA^. Moreover we note that the microcanonical ensemble averages of 
molecular particle numbers also satisfy ETH. Due to the total particle number conversation 
N = 2 (^(a[(^^ai(/?))c, + 2(^«|6{^(^)6L(ij))«j, {b^fi)h(R))a satisfies ETH when (al(R)aL(ij))a 
does ETH. 

In contrast to the chaotic cases, the quantum-mechanical mean values strongly depend 
on the energy eigenstate labeled by a in the regular region. This is clearly shown in Fig. [161 

Next, we investigate whether chaos and ETH appear simultaneously. First, we evaluate 
fluctuations of quantum-mehcanical mean values around a microcanonical-mean value. By 
using the definition of Eq. ([6]), we represent the fluctuation around the ith energy eigenstate 
as 



1 1 i+AAT 2 

\ i'=i-AN 



O 1 i+AAT 2 

\ i'=i-AN 



14 




an isw is«i SOCK) wh mco 



FIG. 13: (color online) The distribution 
in the energy-eigenstates in the initial state 
"casei ". This data is in the chaotic region 
y/Ng/Ja =4. a is the numbering to the 
energy-eigenstate. The energy-eigenstates are 
arranged in ascending order. 
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FIG. 14: (color online) The distribution 
in the energy-eigenstates in the initial state 
"caseS ". This data is in the chaotic region 
^/Ng/Ja = 4. Q is the numbering to the 
energy-eigenstate. The energy-eigenstates are 
arranged in ascending order. 




FIG. 15: (color online) a-dependence of the quantum-mechanical mean values in the chaotic region 
y/Ng/ Ja = 4. a is the numbering to the energy-eigenstate. The energy-eigenstates are arranged in 
ascending order. This figure represents about the atomic-particle numbers in left and right wells. 
The straight line is the average of We set as = 40. 



where we set as AA^ = 10. The maximum values of Eqs. (fT3ll and ( IT4ll are normalized to 
1. If Eqs. ( |T3|) and ( |T^ are small, quantum- mehcanical mean values are almost the same, 
which means the realization of ETH. However, from Eqs. ( IT3l) and (IHj), we can evaluate 
the fluctuation around the iih. energy eigenstate only. Therefore, we define the indicator of 
the verification of ETH as 

^{NaLiaWiY^"^"^ = ^ ^ «L(iJ))™x'(2AAf+l)-A7V] 5 (15) 

where Neth = [{Emax — Emm) I (2AA^ + 1)J . Here, \_x\ is a floor function, which gives the 
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FIG. 16: (color online) a-dependence of the 
quantum-mechanical mean values in the regu- 
lar region \/iV g/ Ja = 0.2. a is the numbering 
to the enregy-eigenstate. In this figure, we fo- 
cus on the atomic-particle numbers in the left 
and right wells. 
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FIG. 17: (color online) The fluctuations 
around the microcanonical mean values of the 
atomic-particle numbers in the left and right 
wells. 



largest integer less than x or equal to x. Here we omit the lower and higher 20% energy 
eigenstates in the same way as Fig. |3l Emax {Emin) is the maximum (minimum) value in 
this energy spectra. We numerically evaluate the verification of ETH by Eq. ( lT5l) . We 
define the indicator A{Ni,L{bR))^''°''"''' for molecular particle number in the left (right) well in 
the same way as Eq. (fTSl) . 

From Fig. [T71 ETH is verified with increasing the atom-molecule internal tunneling. 
In particular, the fluctuations about the microcanonical averages are suppressed when 
VNg/Ja > 1- AiNbLibR))^'"^"^ exhibits exactly the same behavior as Fig. [T71 Comparing 
this result with Fig. [3l we conclude that the onset of chaos and ETH occurs simultaneously. 



V. PARTICLE-LOCALIZED EIGENSTATES 



We show that the atom-molecule internal tunneling change energy eigenstates to be ex- 
ponentially localized. We investigate the case = 40 throughout Sec. |Vl 

We write the interwell relative atomic (molecular) particle number as Za(b) ^ 
[NaL(bL) — NaRibR)) /N . Here we define the distribution function of the interwell relative 
atomic (molecular) particle as 

/a(b) (^a(b)) = ^ \^i{NaL,NaR,NbL,NbR)\^6Nz,(^,yN^L(bL)-NaR{bR)^ (16) 

NaL,NaR,Nl,L,NbR 
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FIG. 18: (color online) The distribu- 

We investigate the 



tion functions of ^^(b)- 
1000th energy eigenstate in the Wigner region 
VNg/Ja = 4. 




FIG. 19: (color online) The logarithm of 
fa^'^^iza) in Fig. [181 We omit the odd num- 
bers of Nza- The straight line is the result of 
the least-square method. 



which is the distribution function of the ith energy eigenstate. 

In Fig. [IHl we show the distribution functions of interwell relative particle numbers in 
the Wigner region y/Ng/Ja = 4. In this figure, we investigate the 1000th energy eigenstate. 
When Nza is an odd number, /* (Nza) is zero because the total particle number is the even 
number. Furthermore, we can see the exponential localization of the distribution functions 
in both (Nza) and (Nzb). In Fig. [191 we plot the logarithm of (A^^^) of the 
1000th energy eigenstate. We omit the points where Nza is an odd number. In this figure, 
by using the least-square method, we can approximate the distribution function by straight 
lines. When we plot the straight lines, we omitted the edge points such as Za = ±1- Because 
of the finite size effect, these points are not approximated well by straight lines. 

Next, we evaluate quantitatively the difference between the distribution function and the 
exponent functions. We evaluate this difference numerically as 



7V-1 



Ri^ Yl H [fa (Nz,)] - (Co + C,Nz, 

NZa=0 



(17) 



where we decide the coefficients Cq and Ci by using the least-squared method. When 

/Af2a=0' 



we calculate J2nz^=o^ omit odd numbers of Nza- In Eq. f lTTll . we use the ith energy 
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FIG. 20: (color online) The difference be- 
tween the energy eigenstates and exponent 
functions defined in Eq. (jl7p and Eq. (jlSp . 



FIG. 21: (color online) The difference be- 
tween the energy eigenstates and exponent 
functions defined in Eq. (jl7p and Eq. (jlSp . 
In this figure, we focus on the Poisson region 
(0 < ^gjJa < 0.1). 



eigenstate. Furthermore, we define the average value over a number of the energy eigenstates. 

Ri 



R 



Ave 



2000 

E 

i=1000 



1001 



(18) 



where we average the RiS from 1000th to 2000th energy eigenstates. When = 40, the 
number of the energy eigenstates is 3311. From 1000th to 2000th energy eigenstates exist 
in the middle of energy spectra. We plot the -Riooo and Ravc in Figs. [SU] and [21] 

Fig. [20] shows that the energy eigenstates behave exponentially when the interwell and 
internal tunneling compete. While -Riooo and Ravs are extremely large in the region < 
VNg/Ja < 0.1 as shown in Fig. [211 -Riooo and Ravc are small in the region VN g ~ Ja as 
shown in Fig. [20] Comparing this behavior with Fig. [3] we conclude that the exponentially 
localization of the energy eigenstates occurs simultaneously with chaos. This means that the 
energy eigenstates behave as exponentially localized functions when thermalization occurs. 
In contrast, the energy eigenstates are quite different from exponential functions in the 
regular region as shown clearly in Fig. [21] 

Next, we compare -Riooo with Ravs- In both Figs. [20] and [211 -Riooo with Ravc do not 
coincide. In particular, the 1000th energy eigenstate does not always localize exponentially 
as shown in Fig. [20] -Riooo fluctuates. The exponential localization clearly appears in the 
averaged value Ravc- This fact indicates that the significant nature of the chaotic eigenstates 
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FIG. 22: (color online) The indicator for the 
localization of energy eigenstates. Lavb is de- 
fined in Eq. ([20]) . 



FIG. 23: (color online) The fluctuations N^ve 
and TAve defined in Eqs. ([MD and ([22]). 



manifests itself by some kind of statistical process. This is the same as the idea of the level 
statistics. 

Next we show that the distribution functions are delocalized when the energy eigenstates 
behave as exponential functions. We define the indicator for localization as 



N 

Nza=-N 



(19) 



When the distribution function of the ith energy eigenstate /* (Nza) is localized in a specified 
Nza only, Lj is unity. When the distribution function of the zth energy eigenstate /* (Nza) 
is spread out over many Nza, Li is close to zero. Furthermore we define the average of LjS 
over many energy eigenstates as 



2000 



-'Ave 



1000 



1001 



(20) 



We investigate L^ve in Fig. From this figure, delocalization occurs in the chaotic region. 
In contrast energy eigenfunctions are localized in the regular region. 

In Figs. [20] and [211 we investigated the exponential localization. By comparing this 
with Fig. [221 we note that the energy eigenfunctions are delocalized when the exponential 
localization occurs. When thermalization occurs, the particles move freely. 

In the last of this section, we investigate the interwell atomic phase fluctuations. The 
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canonical-conhigate variable of the interwell relative particle number Zafb) is the interwell rel- 



ative phase 



33| . When Za(b) behaves as exponential functions, the interwell phase fluctuation 



is increased as shown below. 

We define the fluctuation of the interwell atomic tunneling in the ith energy eigenstate 

as 



T,- = 



N 



\ 



(21) 



where the maximum value is unity. We define the average of TjS over many energy eigenstates 

as 



2600 



Tavc = ^ ] 



601 



2000 



(22) 



Furthermore we define the fluctuation of S^aLa^aR as 



Ni = — 
N 



(23) 



where the maximum value is unity. We define the average of A^jS over many energy eigen- 

SijcltjGS dbS 



N 



Ave 



2600 

E 

601 



2000' 



(24) 



In Fig. [2S1 we investigate T^ve and N^ve- When we replace the creation- annihilation oper- 
ators by c-numbers as dx,(R) ~ NaL{aR)e^^''^'^''^\ {o}l'^r + dj^aj^j /2 becomes 



^/NalNaR COS {OaL-Oa 



R) 



(25) 



In Fig. mi Tavs is increased while N^ve saturates. From Eq. f l25p . this means that the inter- 
well atomic phase fluctuation is increased. The exponential behavior of Za is accompanied 
with the increase of the interwell atomic phase fluctuation. 
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FIG. 24: (color online) The C/a-dependence 
of r/. In this figure, only even spectra are in- 
cluded. The lower and higher 20 % levels are 
not included. We set y/Ng/Ja = 2. 
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FIG. 25: (color online) The fluctuations 
around the microcanonical mean values of the 
atomic-particle numbers in the left and right 
wells. We set VNg/Ja = 2. 



VI. THE INFLUENCE OF THE INTERPARTICLE INTERACTION 

In this section, we instigate the influence of the interatomic interaction on our results. 
Throughout this section, we set \/Ng/Ja = 2. As shown in Fig. [3], chaos occurs at 
y/Ng/Ja = 2 when NUa/Ja = 0. Furthermore, as shown in Fig. [T71 ETH is verified 
simuhaneously. We investigate the influence of the finite Ua on chaos and thermahzation. 

In Fig. [24l we plot rj by increasing NUa/Ja- In this figure we omit the lower and higher 
20% energy eigenstates as in Fig. [31 From this figure, we conclude that the interatomic 
interaction suppresses chaos in our system. When NUa ~ ''7 is near zero. The system is 
chaotic. However, rj is increased by the inter atomic interaction in the large NUa/ Ja region. 
The periodicity of dynamics is recovered. When NUa/ Ja is increased further, rj saturates to 
some constant below unity. This means that the periodicity of dynamics is not completely 
recovered. 

In Fig. [25l we investigate the influence of the Ua on ETH. In Fig. [25], we investigate 
^{N aL{aR))^^°^"'^ ■ As Same as r/, ^{NaL{aR))^^°^°'^ is increased and saturates in the large 
NUa/ Ja region. This indicates that the breakdown of ETH. However, comparing with Fig. 
[T71 the fluctuations /^{NaLiaK))^^"^"^ are small in Figs. [251 This reflects the fact that the 
revival of the periodicity is not complete. We note that ^{Ni^m^^-^)^^"^"'^ exhibits exactly the 
same behavior as A{NaL{aR))^''°'"^'' ■ 
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VII. SUMMARY 



We investigated how the atom-molecule internal tunneling deform the level-spacing dis- 
tribution and dynamics. Increasing the internal tunneling strength, the level-spacing dis- 
tribution changes to Wigner distribution from Poisson distribution. Further increasing the 
internal tunneling, the level spacing distribution becomes Poisson distribution again. The 
atom-molecule internal tunneling changes the dynamics significantly. 

Furthermore, we investigated the thermalization dynamics. This thermalization is ac- 
companied by the manifestation of eigenstate thermalization hypothesis (ETH) 2J]. As 
explained in Sec. IIVB[ ETH explains the principle of equal weight in statistical mechanics 
from the microscopic view point. The foundation of statistical mechanics is based on the 
significant change of energy eigenstates. We showed numerically that the onset of ETH 
occurs simultaneously with that of chaos. 

We also found that the energy eigenstates become exponentially localized as shown in 
Sec. |Vl This exponential localization is analogous to Anderson localization. However, our 
result is significantly different from it in one respect. In the previous studies {2, 3, 5, 34) 35|. 
the localization occurs in the insulator region, where the level statistics obeys the Poisson 
distribution. In contrast, the exponential localization appears in the Wigner region in our 
study. The opposite behavior occurs. 

In our study, interparticle interactions are not necessarily needed for thermalization. 
The atom-molecule internal tunneling induces thermalization in the absence of interparticle 



interactions. This point is different from the previous studies |24j. l28l. |30| . In these studies, 
Wigner distribution is induced in the presence of interparticle interaction. In Sec. I VI I 
only, we include the interparticle interactions in our analysis. As a result, we showed that 
thermalization occurs when NUa, J a and \/Ng competes. It seems that the competition 
between parameters is important for thermalization. 

We note that the multiplicity of deg rees of freedom causes thermalization independently 
of integrability in general systems 23|. When the number of degrees of freedom is large, 
thermalization occurs. This is quite different from ETH. Not the multiplicity of degrees of 
freedom but non- integrability plays an important role in ETH. At least, large system size 
is not needed for ETH. In this paper we show that ETH is verified in the system with only 
two sites. Our study clearly shows that large system size is not important for ETH. 
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Finally, we note that the same result as atom-molecule Bose gases cannot occur in a binary 
Bose gases. By replacing the internal tunneling term in Eq. (|2]) with [b\aL + b^jidR + a^^L + 
a^^i?) ; we obtain the Hamiltonian for a binary BEC mixture. From the same procedure as in 
this paper, we can calculate the level statistics of a binary Bose gases. As a result, it turned 
out that Wigner distribution does not appear. Thermalization is caused by atom-molecule 
internal tunneling. 
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Appendix A: Parameters 



The parameters in the four-mode Hamiltonian ([2]) are defined as follows : 



Ji = - I dT<^lL 



2mi 



2mi 



iL 



2m,- 



Ua=ga / dr\<!),L\^ = Qa / rfr|$ 



■iR\ 



^ = S J rfr|$5L|' + - 2E° = 6 J dr\^,R\^ + E° - 2E°, 



(Al) 

^^R, (A2) 

(A3) 
(A4) 

(A5) 



where i = a (b) represents atomic (molecular) lowest-energy modes, and L (R) expresses the 
left (right) well respectively. 
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